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00 ! We present optical lightcurves of the gravitationally lensed components 

A (= A1+A2+A3) and B of the quadruple quasar RX J0911.4+0551 (z = 2.80). 
The observations were primarily obtained at the Nordic Optical Telescope be- 

^ ; tween 1997 March and 2001 April and consist of 74 /-band data points for each 

component. The data allow the measurement of a time delay of 146 ± 8 days 

O ' (2c) between A and B, with B as the leading component. This value is signifi- 

cantly shorter than that predicted from simple models and indicates a very large 
external shear. Mass models including the main lens galaxy and the surrounding 
Ph' massive cluster of galaxies at z — 0.77, responsible for the external shear, yield 

H = 71 ±4 (random, 2a) ± 8 (systematic) kins -1 Mpc -1 . The systematic model 
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uncertainty is governed by the surface-mass density (convergence) at the location 
of the multiple images. 

Subject headings: cosmology: observations — dark matter — distance scale - 
gravitational lensing — quasars: individual (RX J091 1.4+0551) 



1. Introduction 

The light from a gravitationally lensed, multiply imaged quasar travels along slightly 
different paths towards the observer. This gives rise to a time delay: if the quasar is variable 
the same variation will be seen in the multiple images at different times. The time delay is 
inversely proportional to the Hubble constant, thus offering a means to measure its cosmo- 
logical value in a direct way (Refsdal 1964). Depending on the source and lens redshifts the 
time delay will have a weak dependence on other cosmological parameters. The time delay 
also depends on the gravitational potential and hence the detailed mass distribution towards 
the quasar, thus complicating the measurement of cosmological parameters but providing 
powerful constraints on distant structures such as galaxies, groups, and clusters. 

With the aim of constraining cosmological parameters and the mass distribution in 
galaxies and clusters we have conducted a survey of time delays for multiply imaged (gravi- 
tationally lensed) quasars with the Nordic Optical Telescope (NOT). A description of the pro- 
gram and a measurement of the time delay in our first target, the double quasar B 1600+434 
which is lensed by an edge-on disk galaxy, was presented by Burud et al. (2000). We here 
present the observations, data analysis and time delay of the second survey target, the 
quadruple quasar RX J0911. 4+0551 (z = 2.80) which is lensed by a galaxy and a cluster at 
z = 0.77. The time delay of the third survey target, SBS 1520+530, is presented by Burud 
et al. (2002). 

RX J091 1.4+0551 was discovered by Bade et al. (1997) in the ROSAT All Sky Survey. 
Burud et al. (1998) found optical and near infrared evidence that RX J0911. 4+0551 is a 
quadruply lensed quasar with an unusual image configuration requiring a large external 
shear. The origin of this large external shear was attributed to a possible nearby cluster 
at a photometric redshift z = 0.7 ± 0.1. Kneib, Cohen & Hjorth (2000) confirmed the 



1 Based on observations made with the Nordic Optical Telescope, operated on the island of La Palma 
jointly by Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofisica dc Canarias. 
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existence of a massive cluster and measured the redshift z = 0.769 and velocity dispersion 
a = 836I200 kms -1 of the cluster (to which the main lens galaxy belongs). From Chandra 
X-ray Observatory observations Morgan et al. (2001) found a cluster temperature of 2.31q;s 
keV and a 2-10 keV luminosity of 7.6±g:| x 10 43 erg s _1 . 

The three A (Al, A2, A3) components of RX J0911.4+0551 are very close (0 A i - #A2 = 
0'.'478, 0A2 — &A3 = 0'.'608) and the time delays between them are expected to be short, of the 
order of less than one week. Therefore, we focused on determining the time delay between 
A and B, separated by 3.1 arcsec, expected to be of the order of many months. The images 
are fairly bright with mean /-band magnitudes of 17.2 and 19.2 for A and B respectively. 
Early observations showed that the QSO is strongly variable with a time delay of about 200 
days (Hjorth et al. 2001). We here present an analysis of the full data set. 



2. Lightcurves 

/-band images were obtained at the NOT about once a week between September 1998 
and April 2001. These regular monitoring data were supplemented with a few early data 
points obtained at NOT and MDM between March 1997 and June 1998. The target is 
below the horizon at the NOT in July and August. There were additional gaps in the 
lightcurves due to periods of bad weather and time allocated to Spanish and international 
observing programs. Three different instruments were used: ALFOSC (Andalucfa Faint 
Object Spectrograph and Camera), HiRAC (High Resolution Adaptive Camera) and the 
stand-by camera StanCam, equipped with detectors yielding pixel scales of / /189, 0'.'107 and 
/ .'176 respectively. The /-band was chosen to minimize the sensitivity to lunar phase. One 
data point typically consisted of three exposures of 3-5 minutes each. The seeing varied from 
0'.'6 to l'/5, with 0f9 being the most frequent value. We typically obtained a signal-to-noise 
ratio of 200-300 for the summed A components and 50-100 for the B component. 

The data were reduced and analyzed as described by Burud et al. (2000). Preprocessing 
was done by dedicated pipelines while fringe correction and cosmic-ray removal were per- 
formed manually. Data from different detectors were brought onto the same photometric 
reference system (Burud et al. 1998) via appropriate color terms. Three reference stars with 
known magnitudes were used to calibrate the photometry (see Table 1 for their coordinates 
and magnitudes). The photometry of the quasar images was performed by applying the 
MCS deconvolution algorithm (Magain et al. 1998). This algorithm has already been used 
to analyze the data of several blended lensed quasar images (e.g., Burud et al. (1998), Burud 
et al. (2000)). Its main advantage is its ability to use all, even rather poor, data, irrespective 
of image quality and lunar phase. The final deconvolved image is produced by simultane- 



-4- 



ously deconvolving the individual frames of the same object from all epochs. The positions 
of the quasar images and the shape of the lensing galaxy are the same for all the frames 
and are therefore constrained using the total S/N of the data-set. The intensity of the point 
sources are allowed to vary from image to image, hence producing the lightcurves. Images 
of RX J091 1.4+0551 and the resulting lightcurves are presented in Figs. 1 and 2. 

3. Time delay 

The J-band lightcurves (Fig. 2) contain 74 data points for each component. As predicted 
by theory, the lightcurves show that B is the leading component. A pronounced V-shaped 
feature at JD 2451300 (May/June 1999) is seen in the A component followed by several 
decreases and upturns. These are preceeded by similar features in the B component which 
allows the determination of a rough time delay of about 150 days. The data points are 
available upon request. 

A quantitative analysis of the light curves was performed using the four methods de- 
scribed by Burud et al. (2000). The SOLA method (Pijpers 1997) does not provide a definite 
time delay but is consistent with the results of the minimum dispersion method (Pelt et al. 
1996) and the two novel methods introduced by Burud et al. (2001). The time-delay esti- 
mates and magnitude offsets obtained from the three different methods are consistent with 
each other and presented in Table 2. 

It is readily apparent that no simple time translation will turn the A curve into the B 
curve. Thus, 'external' variations in the time-delay corrected flux ratios must be present. 
This is confirmed by our analysis of the individual Al, A2 and A3 lightcurves which exhibit 
similar overall trends but different detailed shapes. The time delays determined from the 
individual components are consistent with the expectation of small time delays between the 
A components and with the average B-A time delay (see Table 3). 

The 'external' variations, due e.g. to microlensing (see Burud et al. (2000)) are best taken 
into account by the methods introduced by Burud et al. (2001). In fact, the cause of the 
failure of the SOLA method lies in the nature of the external variations in RX J0911. 4+0551, 
which are of bigher order than linear in time over the entire time series. This violates an 
assumption under which that method is derived. The best method for dealing with such high- 
order variations is the iterative method. We therefore adopt the time delay as determined 
by the iterative method (see Fig. 3) and conclude that the time delay between A and B in 
RX J0911.4+0551 is 146 ± 8 days (2a). With a 2a error of 5 % this is among the most 
precise time delays determined for any lens system. 
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4. Discussion 

In modeling the system we used a cosmology with Q m = 0.3, Q\ = 0.7, and Hq = 
lOO/i kms -1 Mpc -1 . Adopting an open Universe with VL m = 0.3, £l\ = increases the model 
time delay (and derived H ) by 4.1 percent, whereas an Einstein-de Sitter Universe with 
VL m — 1.0, Q,\ — leads to a decrease of 10.0 percent. 

The 'yardstick' potential model of Schechter (2000) involves a two-dimensional potential 
given by ip ~ 6 1+a with a = —0.8 and an external shear with k = 7 = 0.307. In the adopted 
cosmology, this model predicts a time delay of 1 1 1/i 1 with an estimated uncertainty of 
±22/i~ 1 days due to the uncertainty in a from the scatter in the observed shapes of nearby 
elliptical galaxies. In the model presented by Kneib, Cohen & Hjorth (2000) the predicted 
time delay is 112.5 ± 17.5/i _1 , consistent with the 'yardstick' prediction of Schechter (2000), 
and with a slightly smaller uncertainty. 

We independently modeled the system as described by Kneib, Cohen & Hjorth (2000). 
The model includes the main lensing galaxy, the cluster of galaxies, and individual galaxies 
in the cluster. We used the measured velocity dispersion of Kneib, Cohen & Hjorth (2000) 
to constrain the cluster mass and adopted more realistic contraints on the ellipticity of the 
cluster mass distribution and a tighter mass-luminosity relation for the cluster galaxies. The 
main model uncertainty concerns the value of the cluster convergence k at the location of the 
multiple images. A large k gives rise to a small predicted time delay and derived H . The 
value of k is governed by the mass of the cluster, its mass profile, and the effects of possible 
substructure in the cluster. We find that 0.20 < k < 0.28 is required for a good fit. This 
range in k translates into different values for the mass and velocity dispersion of the cluster, 
depending on the mass profile used and the ellipticity of the cluster mass distribution. In 
addition to the smooth cluster convergence there is a contribution to the total convergence 
of about 0.06 from the individual galaxies in the cluster (not including the main lens). 

The refined model prediction of the (flux-weighted or straight) mean B-A time delay 
is 104 ± 1 1 (the individual time delays between the A components are less than 1.5 
h^ 1 days, generally in the sequence A2, Al, A3). Using the measured B-A time delay of 
146 ± 8 days (2a) , the resulting value of the Hubble constant is H = 71 ± 4 (random, 2a) ± 
8 (systematic) kms^ 1 Mpc -1 . 

The variability of RX J091 1.4+0551 is sufficiently strong and erratic that the prospects 
for refining the time delay to better than ±5 percent are very good. Measurements of the 
time delays between Al, A2, and A3 also appear within reach from intensive optical or X- 
ray monitoring (Chartas et al. 2001). Moreover, further mapping of the cluster potential 
towards RX J091 1.4+0551 from Chandra X-ray Observatory, XMM-Newton, HST, and 
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VLT observations will help bring down the systematic model uncertainties by determining 
the cluster convergence at the location of the QSO images. Finally, the fairly high redshifts 
of the lens and source result in a sensitivity of the order of 10 percent to the adopted world 
model. With a smaller systematic uncertainty in the model and independent constraints on 
Hq the system may be used to contrain Q m and A. Thus, RX J0911. 4+0551 appears as 
one of the most useful individual lens systems for cosmological parameter determination and 
studies of the mass distribution in galaxies and clusters. 

We thank the NOT Director Vilppu Piirola for granting us observing time for this 
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who have contributed to this project by performing the scheduled observations. This project 
was conceived in 1997 while JH, AOJ, JPK, and JP were visiting scientists at the Center 
for Advanced Study in Oslo. We thank Rolf Stabell and Sjur Refsdal for inviting us there 
and for their kind hospitality. JH appreciates the hospitality of the OMP where some of this 
work was conducted. The project was supported by the Danish Natural Science Research 
Council (SNF). IB was supported in part by contract ARC94/99-178 "Action de Recherche 
Concertee de la Communaute Frangaise (Belgium)" and Pole d'Attraction Interuniversitaire, 
P4/05 (SSTC, Belgium). PLS acknowledges the support of US NSF grant AST-9616866. 
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acknowledges financial support by the Theoretical Astrophysics Center (TAC), a collabo- 
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Table 1. Coordinates and /-band magnitudes for three reference stars in the field 

surrounding RX J091 1.4+0551. 





RA(J2000) 


Dec(J2000) 


I 

(mag) 


Rl 


09:11:21.16 


+05:50:43.2 


17.34+0.01 


R2 


09:11:26.54 


+05:51:43.7 


16.30+0.02 


R3 


09:11:29.08 


+05:49:30.0 


17.07+0.01 



Table 2. Estimated time delays and /-band magnitude differences for RX J0911. 4+0551 
calculated with three different methods (Pelt et al. 1996; Burud et al. 2001). 





At 


Am 




(days) 


(mag) 


Minimum dispersion 


153+3 




x 2 fit 


147+4 


1.88 + 0.01 


Iterative fit 


146+4 


1.95-2.05 



Note. - The quoted uncertainties in the 
time delay estimates are la errors. 
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Table 3. Estimated time delays for RX J0911. 4+0551 calculated using the iterative 

method (Burud et al. 2001). 



At 

(days) 

B-A 146±4 

B-Al 143±6 

B-A2 149±8 

B-A3 154±16 



Note. - The 
quoted uncertain- 
ties in the time 
delay estimates 
are la errors. 
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Fig. I. — Stacked /-band images of RX J0911. 4+0551 from a total of ~ 7 hours of exposure. 
The left panel shows the combined image (seeing = 1'.'20). The right panel shows the decon- 
volved image (FWHM = 0'.'28). The main lensing galaxy can be seen close to the three A 
components. North is up and East is to the left. Each frame measures 12 x 12 arcsec. 
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Fig. 2. — /-band lightcurves for RX J0911. 4+0551 A,B and two of the stars in the field 
around RX J0911. 4+0551 that were not used to construct the PSF. The magnitudes are 
calibrated using reference stars in the field (Table 1). The plotted B lightcurve is offset by 
— 1.5 mag and SI is offset by —0.2 mag. The error bars represent the combined photon noise 
and PSF errors, determined as described by Burud et al. (2000). 
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Fig. 3. — Combined lightcurve from both components of RX J0911. 4+0551. The curve of 
the A component is shifted forward in time by 146 days. The B curve has been offset by 
values between —1.95 mag and —2.05 mag. 



